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Abstract 

We study monomial ideals, always locally given by a monomial, like a reasonable first step 
to estimate in general the number of monoidal transformations of Villamayor's algorithm of 
resolution of singularities. The resolution of a monomial ideal < X^ 1 ■ . . . ■ > is interesting 
due to its identification with the particular toric problem < Z c — X" 1 • . . . • X„ n >. 

In the special case, when all the exponents are greater than or equal to the critical 
value c, we construct the largest branch of the resolution tree which provides an upper bound 
involving partial sums of Catalan numbers. This case will be called "minimal codimensional 
case". Partial sums of Catalan numbers (starting 1,2,5,...) are 1,3,8,22,... These partial 
sums are well known in Combinatorics and count the number of paths starting from the root in 
all ordered trees with n + 1 edges. Catalan numbers appear in many combinatorial problems, 
counting the number of ways to insert n pairs of parenthesis in a word of n + 1 letters, plane 
trees with n + 1 vertices, . . ., etc. 

The non minimal case, when there exists some exponent a,: smaller than c, will be called 
"case of higher co dimension" . In this case, still unresolved, we give an example to state the 
foremost troubles. 

Computation of examples has been helpful in both cases to study the behaviour of the 
resolution invariant. Computations have been made in Singular (see [8]) using the desing package 
by G. Bodnar and J. Schicho, see [1]. 

1 Introduction 

The existence of resolution of singularities in arbitrary dimension over a field of characteristic zero 
was solved by Hironaka in his famous paper [5j. Later on, different constructive proofs have been 
given, among others, by Villamayor [13], Bierstone-Milman [T], Encinas-Villamayor [BJ, Encinas- 
Hauser [5] and Wodarczyk [14]. 

This paper is devoted to study the complexity of Villamayor's algorithm of resolution of singu- 
larities. This algorithm appears originally in [T3] and we will use the presentation given in [BJ. In 
this paper, the authors introduce a class of objects called basic objects B = (W, (J, c), E) where W 
is a regular ambient space over a field k of characteristic zero, J C Ow is a sheaf of ideals, c is an 
integer and E is a set of smooth hypersurfaces in W having only normal crossings. That is, they 
consider the ideal J together with a positive integer c, or critical value defining the singular locus 
Sing{J, c) = {£ € W\ ord^(J) > c}, where ord^(J) is the order of J in a point £. 

Let W W be the monoidal transformation with center Z C Sing(J,c), tt^ 1 {Z) = Y' is the 
exceptional divisor. Let £ be the generic point of Z, ord^(J) = 8, the total transform of J in W 
satisfies JOw = I{Y') e ■ J r where J Y is the weak transform of J, (see [BJ for details). 
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A transformation of a basic object (W,(J,c),E) <— (W , (J' , c), E') is defined by a monoidal 
transformation W W and defining J' = I(Y') e ~ c ■ J Y , the controlled transform of J. 
A sequence of transformations of basic objects 

(W,(J,c),E) «- (WW,(jW,c),EW) - ••• - (^)(J< N >c),#)) (1) 

is a resolution of (W, (J, c), J5) if Sing(j( N \c) = 0. 

Remark 1.1. Superscripts ( fe ) in basic objects will denote the fc-stage of the resolution process. 

(k) 

Subscripts j will always denote the dimension of the ambient space W> . 

Villamayor's algorithm provides a log-resolution in characteristic zero. A log-resolution of J is 
a sequence of monoidal transformations at regular centers as ([T]) such that each center has normal 
crossings with the exceptional divisors E^' , and the total transform of J in W^ N ' is of the form 

JO w{N) = I{Hx)^ • . . . • I(H N ) b » 

with b t e N for all 1 < i < N and E^ = {Hi, H N }. 

In [5J it is shown that algorithmic principalization of ideals reduces to algorithmic resolution of 
basic objects. That is, starting with c~max ord(J), the maximal order of J, we obtain a resolution 
of (W, (J,c),E) as ([TJ. At this step max or d(jW)= c W< c. If C W> 1, we continue resolving 
(W<- N \ (JW,cW), and so on, until have max ord(jW)= c W = i. Finally, a resolution of 

(W^\ (J^\ 1), fiW) provides a log-resolution of jW, and therefore a log-resolution of J. 

In [5] it is also shown that algorithmic principalization of ideals leads to embedded desingulari- 
zation of varieties. That is, given a closed subscheme X C W, the algorithmic principalization of 
the ideal I{X) provides an embedded desingularization of X. See also [7] for more details. 

A key point in the definition of the algorithm is to use induction on the dimension of the ambient 
space W to define an upper-semi-continuous function t. The set of points where this function attains 
its maximal value, Max t, is a regular closed set, and defines a regular center for the next monoidal 
transformation. 

A resolution of the basic object (W, (J,c),E) is achieved by a sequence of monoidal transfor- 
mations as in (QJ, with centers Max t^ for < k < N — 1. That is, the sequence of monoidal 
transformations is defined by taking successively the center defined by the upper-semi-continuous 
function. The algorithm stops at some stage because the maximal value of the function t drops after 
monoidal transformations, that is, max > max t y ' > . . . > max 

This function t will be the resolution invariant. We shall work with the invariant defined in [B], 
using the language of mobiles developed in [5]. We remind briefly the main notions. 

Let J C Ow be an ideal defining a singular algebraic set X C W. The ideal J factors into 
J = M ■ /, with M the ideal defining a normal crossing divisors, and I some ideal still unresolved. 

By induction on the dimension of W , we will have this decomposition at every dimension from n 
to 1, that is Ji = Mi -Ii, for n > i > 1, are defined in local flags W n 3 W n -i D ■ ■ ■ D Wi D ■ ■ ■ D Wi, 
where each Ji,Mi,Ii £ Owi are in dimension i. There is a critical value Ci+\ at each dimension i, 
(c n +i = c), see [5] for details. All the basic objects (Wi, (Ji, Cj + i), Ei), for n > i > 1, will be resolved 
during the process of the algorithm. 

Let E be the exceptional divisor of previous monoidal transformations, and consider E = \Jf =1 Ei 
where Ei applies to dimension i. Obviously, we start with E = 0. 

For any point £ £ Sing(J 7 c), the function t will have n coordinates, with lexicographical order, 
and it will be one of the following three types: 



(a) = (i n (£),i n _i(f),...,t n _ r (f), oo, oo,...,oo) 

(b) t(0 = (*n(0.*n-l(0.---.*n-r(0.r(0.OO,...,Oo) With ti = 
(C) *(0 = (*»(0.*n-l(0,...,*»-r(0. .*l(0) 



^+1 



(2) 
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where 6i = ord^(Ii) , mi is the number of exceptional divisors in E i: and T is the resolution function 
corresponding to the so-called monomial case, following the notation of [6], pages 165 — 166. We 
will recall the definition of T in equation ([JJ . 

For simplicity, let assume that we start with a polynomial ring, W = Spec{k[X\, . . . , X n \). In 
Ow = k[X\, . . . , X n ] the ideal J is locally given by a monomial with respect to a regular system of 
parameters 



J =< I" 1 • . . . • A£" >C O w with a t e N, for i = 1, . 



, n. 



Note that, in this situation, the center of the next monoidal transformation is combinatorial, it 
is a linear combination of X\, . . . , A„. And this is also true after monoidal transformations, since 
Villamayor's algorithm applied to a monomial ideal provides always combinatorial centers, and after 
a monoidal transformation in a combinatorial center we obtain again a monomial ideal. 

So, at any stage of the resolution process, W — UjC/j, where Ui = AJJ. Thereafter, we shall work 
locally, so we will assume that W is an affine space. 

To resolve the toric hypersurface {/ = 0} = {Z c — A" 1 ■ . . . ■ X® n = 0} we note that its singular 
locus Sing(< / >, c) is always included in {Z — 0}, so we argue by induction on the dimension and 
reduce to the case where the corresponding ideal J is of the form 

n 

J =< I" 1 • . . . • X^ >C Ow with 1 < ai < a 2 < . . . < a n , a l = d, d > c, (3) 

i=l 

where c is the critical value. If a, = for some i, then we may assume dim(W) < n. 

After a monoidal transformation, we always consider the controlled transform of J with respect 
to c, J' = I(Y')~ C ■ J* where J* is the total transform of J and Y' denotes the new exceptional 
divisor. For the toric problem J =< Z c — A" 1 • . . . • X°" >, taking the origin as center of the next 
monoidal transformation, at the i-th chart: 

J* =< Z c ■ X c t - XI 1 ■ ■ ■ Xf ■ ■ ■ X% n >=< X? ■ (Z c - X? 1 • • • Xf~ c ■■■X^)>, 

and we can only factorize c times the exceptional divisor. 

Remark 1.2. We will denote as i-th chart the chart where we divide by Xi. When the center of the 
monoidal transformation is the origin, this monoidal transformation is expressed: 

k[Z,X u ...,X n ] -» /■• /.. X V Xi tiM 2L 



Z 



n ' Xi ' Xi ' ' ' ■ ' Xi ' Xi ' ■ • • ' Xi J 

z 

Xi 

Xi — > Xi 

X 3 -> ^ for j ^ i 

where k[Z,X u ■ , ^] = 

k[-£-, . . . , — j^, Xi, —^r, ■ ■ ■ , -j^]- For simplicity, we will denote each jf- again as Xj, and -4- 
as Z. 

So we will apply the resolution algorithm to the basic object (W, (J, c), 0) for J =< Xj" 1 • . . . ■ 
1°™ >, which is already a monomial ideal, but it is not necessarily supported on the exceptional 
divisors. 



2 Monomial case (exceptional monomial) 

The monomial case is a special case in which J is a "monomial ideal" given locally by a monomial 
that can be expressed in terms of the exceptional divisors. This case arises after several monoidal 
transformations. 
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This means that we have a basic object (W, (J, c), E) where J is locally defined by one monomial 
supported on the hypersurfaces in E. In this case, the ideal J factors into J = M ■ I with J = M 
and 1=1. We can also call it exceptional monomial. 

Theorem 2.1. Let J C Oyy be a monomial ideal as in equation f3|). Let E — {Hi, . . . , H n }, with 
Hi = V{Xi), be a normal crossing divisor. 

Then an upper bound for the number of monoidal transformations to resolve (W, (J, c), E) is given 
by 

d — c + gcd(ai, . . . ,a n , c) 
gcd{a 1 ,. . . ,a n ,c) 

Proof. We may assume that the greatest common divisor of the exponents and the critical value 
c is equal to 1, because both the simplified problem and the original problem have the same singular 
locus. That is, if gcd(a 1 , . . . ,a n ,c) = k then d = k ■ rfi, c — k ■ C\, ai = k ■ for all 1 < i < n and 
gcd(b\, . . . ,6„,ci) = 1. The ideal J can be written as J = {J\) k where J\ —< X\ x ■ .. . ■ > 
therefore 

Sing(J,c) = {£ € X\ ord i ((Ji) k ) > k ■ a} = {£ e X\ ord ( (Ji) > ci} = Sing(Ji, Cl ), 

where X is the algebraic set defined by J. 

For a point £ e A" - (-T^),^),^)) where 

r x (0 = min{p | 3 i u ...,i p ,a il (0 + --- + a ip (£) >c, H n n---DH lp }, 

T 2 (0 = max j a ' l(0+ ' e ' +a ' p(0 | p=rx(5), a n (?)+-+a lp (C)>c, Sefii.n-nHi,,}, (4) 

r 3 (0 = nua{(ii, . . . ,i p , 0, . . . ,0) e Z n I r 2 (0= a,l(O+ ; +aip< ° , ee^n-n^} 



with lexicographical order in Z™ . 

The center Z of the next monoidal transformation is given by the set of points where T attains 
its maximal value. It is easy to see that Z = ri^ =n _^ r _ 1 - ) H i . 

So at the j-th chart, the exponent of Xj after the monoidal transformation is {Y27=n— r+i a *) — c 
and 



E 



a i ) — c < mm dj = a n - r +i 

n— r+l<z<n 



\z— n — r+1 / 

because X^™=n-r+2 ai < c by construction of the center Z. 

This shows that the order of the ideal drops after each monoidal transformation by at least one, 
so in the worst case, we need d — (c — 1) monoidal transformations to obtain an order lower than 
c. □ 

Remark 2.2. Note that it is necessary to consider the monomial case. On one hand, this case may 
appear in dimension n, and also in lower dimensions, n — 1, . . . , 1, when we resolve any basic object 
(W, (J,c),E) (where J is any ideal). So we need to resolve the monomial case in order to obtain a 
resolution of the original basic object (W, ( J, c), E). 

On the other hand, the algorithm of resolution leads to the monomial case, since given any ideal 
J, the algorithm provides a log-resolution of J. And it is necessary to continue to a resolution within 
the monomial case. 

Remark 2.3. The bound in theorem 12. H is reached only for the following values of c: 

1, a n + . . . + Oj + 1 for n > j > 2, d. 
For these values of c, the order of the ideal drops after each monoidal transformation exactly by one: 
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• If c = 1, the monoidal transformation is an isomorphism. The exponent of X n after the 
monoidal transformation is a n — 1. 

• If c = a n + . . . + dj + 1, for n > j > 2, the center of the monoidal transformation is Z = 
n" = j_ 1 Hi. At the i-th chart, for n > I > j — 1, the exponent of A; after the monoidal 
transformation is (Y%=j-i <k) —c = (Y%=j-i a i) ~ a i) ~ 1 = a j-i ~ l - 

In particular, at the (j — l)-th chart, the exponent of after the monoidal transformation 
has droped exactly by one. 

• If c = d, we finish after only one monoidal transformation. 

Remark 2.4. If gcd(a\, . . . , a n , c) = k > 1, then the bound of the theorem 12. H is (e? — c + 

As (<i — c + < d — c + 1, we can use in practice the bound for the case gcd(a\, . . . , a n , c) = 1. 

3 Case of one monomial 

To construct an upper bound for the number of monoidal transformations needed to resolve the basic 
object [W, (J, c),E = 0), where J is locally defined by a unique monomial, we estimate the number 
of monoidal transformations needed to obtain (W , ( J' , c), E'), a transformation of the original basic 
object, with J' = M' (an exceptional monomial), and then apply theorem 12.11 In order to use 
theorem 12. 1[ we need an estimation of the order of M' . This estimation will be valid at any stage 
of the resolution process. 

Lemma 3.1. Let (W, (J, c), 0) be a basic object where J is a monomial ideal as in equation f3J). Let 
J = M ■ I be the factorization of J , where M = 1, because of E = 0, and J = I . After N monoidal 
transformations we have {W^ N \{J^ N \c),E^). Let^EW^ be a point. Then 

ord^(M^) < (2 N - l)(d-c) 

where ord^(M^ N ^) denotes the order at £ of M^ N \ the (exceptional) monomial part of J^ N \ 

Proof. It follows by induction on Af: 

• that if N = 1, ord^M^) = d - c. 

At the beginning, the first center defined by this algorithm is always the origin, so at the i-th 
chart: 

J« = AfW . /(i) =< X f- C > ■ < X" 1 ' .?. > 
with E^ = {Hi}, where H, = V(Xi). 

• We assume that the result holds for N = m — 1 . 

j(m-l) = M (m-1) . ^m-l) =< ^ . ..^6. > . < Jf *l.+i . . . j*, > 

with X)i=i ^« = inductive hypothesis, after m— 1 monoidal transformations, the maximal 

order d' of the (exceptional) monomial part Af(™ 1-1 ) satisfies 

tf < (2™- 1 - l)(d-c). 

For N = m, there are two possibilities: 

1. If max ord(L( m ~ lS> ) = X)j= s +i a h — c then the center of the next monoidal transformation 
contains only variables appearing in /t" 1 " 1 ). 
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2. If max ord(l( m ^ ) = Y^j= s +i a; H < c then the center of the next monoidal transformation 
contains variables appearing in I^ rn ~ 1 "> and also variables appearing in M^ m_1 ^. 

Case 1 : 

If the center of the monoidal transformation is as small as possible, that is Z = r\^ =s+1 V{X ij ), 
at the i/-th chart, 

j(m) =M (m) . 7 M =< Vbi. . . X b. . X d "^ =lCS " C >-<X ais+1 >. 

fc l "l ''s-l-l »n 

The exponent of X in d — a ij ~ c = Y^j= s +i ai i ~ c 1S as ki§ as possible, so this is 

the worst case, because the increase in the order of the exceptional monomial part after 
the monoidal transformation will be greater than that for another centers. 
The highest order of is 

s s s 

hi + d — a ij — c = d' + d — c — ^ a ij < d' + d — c, 
*=i i=i i=i 

so by inductive hypothesis 

d' + d - c < (2™- 1 - \){d-c)+d-c = 2 m -\d - c) < (2 m - l)(d - c). 
Case 2 : - At the ij-th chart, for 1 < j < s 

j(m) = M (m) . j{m) =< j^i ^6. . V □ > . < v"^ 1 . . . > . 

- At the ii-th chart, for s + 1 < I < n 

j(m) = M (m) . jim) =< j^i . . . ^6. . yA . ya« B > _ 

As above, if we are in the worst case, when the center of the monoidal transformation is 
as small as possible, that is, the center is a point, 

s 

□ =A= d' + d - y~] - c . 

Therefore in both cases the highest order of satisfies 

< 2d' + d - c < 2(2 m ~ 1 - l)(d - c) + d - c = (2 m - l)(d - c) . 

□ 

Remark 3.2. Due to its general character, this bound is large and far from being optimal. 

Remark 3.3. The ideals Mj are supported on a normal crossing divisors Di. Recall that their 
transformations after monoidal transformations, in the neighbourhood of a point £ G Wi, are 

D* + (6i - Ci+i) • Y' if (t' n (t'),...,t' i+1 (i')=(tn(t),-,U +1 (S)) 



in other case ,»>*>!, 



(r>^=L>*+(e„-c)-y' always) 

where D* denotes the pull-back of by the monoidal transformation ir, Y' denotes the new ex- 
ceptional divisor, the point £' <E W[ satisfies 7r(£') — di ~ ord^(Ii) and c i+ i is the corresponding 
critical value. 
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In what follows we will define the ideals Ji-\, n > i > 1. We need some auxiliary definitions: 
the companion ideals Pi and the composition ideals Ki, see [5] for details. 
We construct the companion ideals to ensure that Sing(Pi,8i) C Sing(Ji, Cj+i) , 

f if #i > Ci+i 

[ h + Mi 1+1 1 if o < e t < Cl+ i 

where £ £ A£ is a point, #i = ord^(Ii) and Cj+i is the corresponding critical value. 

Let Ji — Mi ■ Ii be the factorization of an ideal Ji in Wj , where Mi , Ii are ideals in Wi in the 
neighborhood of a point £ 6 AJJ. Let £^ be a normal crossing divisor in A£. 

The composition ideal i<Q in Wi of the product Ji = Mi ■ Jj, with respect to a control Cj + i, is 

^ _ J Pi -/wi^nwi) if Ji^i, 

Ai ~ \ l if h = l. w 

The critical value for the following step of induction on the dimension is Ci = ord^(Ki). 

The construction of the composition ideal Ki ensures normal crossing with the exceptional divisor 

Ei. 

We say that an ideal K is bold regular if K =< X a >, K e a G N. 

Finally, construct the junior ideal Jj_i 

j _ f Coef fv(Ki) if is not bold regular or 1 , . 

4-1 [1 otherwise 

where V is a hypersurface of maximal contact in Wi (see [5] page 830) and Coef fv{Ki) is the 
coefficient ideal of Ki in V (see [5] page 829). The junior ideal J;_i is an ideal in this suitable 
hypersurface V of dimension i — 1 . 

If ^ > 1 we arc in the first case of equation ([5]), ^— 1 = ^— ^ = . . . = ^— = 1 and tj-i = . . . = 
ti = oo for n — 1 > j > 1, because D n _i — . . . = D\ = and P; = Ij, and hence J,_i is always 
given by a unique monomial. 

Remark 3.4. For an ideal J =< I" 1 • . . . • X°™ > as in equation ([3]), if we assume a n > a n -i > ■ ■ ■ > 
a\ > c, then at every stage ^- > 1 , so we are always in the previous situation. The singular locus 
of ( J, c) is always a union of hypersurfaces U[ =1 {Xj = 0}, l<r<n, and the center of the next 
monoidal transformation will be the intersection of some of these hypersurfaces. So we will call this 
case the minimal codimensional case. 

Remark 3.5. If there exists some a, < c, at a certain stage of the resolution process it may occur 
2°- < 1. Then we are in the second case of equation ([5]), the (exceptional) monomial part M n can 
appear in some Jj for n — 1 > j > 1, and ^— -j- can be much greater than 1, what increase the number 
of monoidal transformations. Now its singular locus is a union of intersections of hypersurfaces of 
the type U; . ({X^ = 0} fl . . . (1 {Xi t — 0}). This is the higher codimensional case. 



4 Bound in the minimal codimensional case 

Remark 4.1. From now on, we always look to the points where the function t, defined in ([2|), is 
maximal. So the following results concerning the behaviour of the function t always affect the points 
where it reaches its maximal value. 
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Proposition 4.2. Let (W,(J,c),E) be a basic object where J is a monomial ideal as in equation 
(0), with a,i > c for all 1 < i < n. We can factor J = J n = M n ■ I n , and after r — 1 monoidal 
transformations, J n r : ' = M n r 1 ' -/i r Let £ £ Vl/( r_1 ) &e a point where ord^(I n r ) —0 n . After 
each monoidal transformation n, the resolution function in a neighbourhood of £' , where 7r(£') = £ 
and ord^ (In ) = 9' n < 8 n , is of the form 

,[1,0],..., [1,0] J for some 1 < s < n - 1. 
Proof. After monoidal transformations, 

jW = mW./« =<^...^>.<^...^- > 

with d — X)j=i a «j = Sj= s +i a «j — c th° n 7 = ^n*^ an d the (exceptional) monomial part does 
not appear in Jj ; for all n > £ > 1. 

We have 0' n ^ 9 n , then .E^ = Y"' + |-E| Y and m„ = s, we count all the exceptional divisors of 
the previous steps and the new one. There are no exceptional divisors in lower dimension because 
= (Y 1 + \E\ Y ) - E n r) = and, in a similar way, we obtain £, (r) = for all n - 1 > I > 1. 

The normal crossing divisors D^' = for all n — 1 > i > 1, so the corresponding ideals 
M n r \ = ... = M[ r) = 1. In particular, M^2i = 1, hence 

c'„ = ord e ,(^ r) ) = ord e (Coeff(KP)) = ord e (j£\) = ord?^) = B' n _ x 

with £' € such that 7r(£') = £, because ord(Coef f(K)) = ord(K) when if is a monomial ideal, 
therefore ^T 1 = 1. By the same argument we obtain 9 ,"~ 2 = . . . = = 1. □ 

Remark 4.3. After each monoidal transformation, the exceptional divisors at each dimension are: 

E ' = { E J if (*»K').-.*5+i«'))=(*n(0,-,ti +I («)) «nd«;=«j 

■? 1 (Y' + (BiU...UE„) Y )-(E^ ^j'+l) in other case 

(E' n =E r if 6^=0,, or B;,=Y' + (_EiU...U£;„) y otherwise) 

for n. > j > 1, where denotes the strict transform of Ej by the monoidal transformation 7r, 
y' denotes the new exceptional divisor, the point £' £ W[ satisfies 7r(£') = £, = ord^'(I'j) and 
0j = ordt(Ij). We denote |£J| = £a U . . . U S n . 

Hence, after the first monoidal transformation, since 9' n < 9 n we have E' n = Y' and E' n _ 1 = 
■ ■ ■ = E[ =0. After the second monoidal transformation, at the chart where = 0' n we obtain 
E% = (E' n ) Y = 0, E'^_ x = Y", and E%_ 2 = ■ ■ ■ = E'{ = and so on. We call this phenomena 
propagation because every exceptional divisor appears in the resolution function t firstly in dimension 
n, then in dimension n — 1, n — 2, and so on. 

Definition 4.4. We will call propagation, p(i,j), for 1 < i < j — 1, 1 < j < n, to the number 
of monoidal transformations needed to eliminate i exceptional divisors in dimension j, when we 
remain constant (t n , t n —i, . . . , £j+i) and 8j, and there are no exceptional divisors in lower dimensions 
j — 1, . . . , 1. That is, passing from the stage 

([6> n ,m„], . . . , [6 j+1 ,m j+1 ], [0j,i], [1,0], . . . , [1,0]) 

to the stage 

i 

([6 n ,m n ], . . . , [6 j+1 ,m j+1 ], [9j,0], [1,0], . . . , [1, 0], bo, , od). 



-,s 



8 



Lemma 4.5. Propagation Lemma Let (W,(J,c),E) be a basic object where J is a monomial 
ideal as in equation with a\ > c for all 1 <l < n. Letp(i,j) be the propagation of i exceptional 
divisors in dimension j in the resolution process of (W, (J, c), E). 
Then, for all 1 < j < n, 



Proof. • If there are i exceptional divisors in dimension i, Ki + \ is bold regular, ti = oo then 
p(i, i) = 0. We can not propagate these i exceptional divisors at this stage of the resolution 
process. 

If there are s exceptional divisors at this step of the resolution process, then there are n — s 
variables in /„. On the other hand, from dimension n to dimension i + 1 there are s — i 
exceptional divisors. 

When we construct J n -i, ■ • ■ , Ji+i, add to the corresponding composition ideal Kj the variables 
in 1^ (Ej n Wj), so in these dimensions there are (n — s) + (s — i) = n — i variables. 

When we make induction on the dimension, at each step lose one variable, so in n— i — 1 steps 
obtain that -fQ+i, that corresponds to the n — (n — i — 1) = i + 1 position, is bold regular. And 
the variables appearing in these i exceptional divisors do not appear in the center of the next 
monoidal transformation. 

• By induction on the dimension: 



- If j = 1, p(l, 1) = by the previous argument. 

- If j = 2, p(l, 2) = 1 because when we propagate 1 excepcional divisor from dimension 2 
to dimension 1, K' 2 is bold regular. 




(8) 



([0 Bj m n ],..., [0 2) 1], [1,0]) 
([6> n ,m„], . . . , [0 2 ,O],oo) 



Then p(l, 2) = 1 = 1 + = 1 + p(l, 1). 
- We assume that the result holds for j < 8 — 1. For j = s: 



([0 n , m„], . . . , [0 s+1 ,m s+ i], [0 s ,i], [1, 0], . . . , [1, 0]) 
I 

([6 n ,m n ], [0 s+ i,m s+1 ], [0„t - 1], [1, 1], [1, 0], . . . , [1,0]) 



I 1 



([0„, m„], . . . , [^+1,7715+1], [0 s ,i - 1], [1,0], ... , [1,0], oo) 



I 

([0 n ,m„], . . . , [0 s +i,m s +i], [0 S , i - 2], [1, 2], [1, 0], . . . , [1,0]) 



We want ([0„,to„] . . . [0 s +i, m s+ i]) and S remain constant. So after the first monoidal 
transformation look to some suitable chart where m s = i drops. As m s drops then 
to s _i = i — (i — 1) = 1 and propagate this exceptional divisor in dimension s — 1, making 
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p(l, s — 1) monoidal transformations. Otherwise, to keep ([0 n ,m n ] . . . [0 s +i, m s +i\) and 
9 S constant, the only possibility is to look to a suitable chart where m s drops from i — 1 
to i — 2. But in this case this would provide the same resolution function that appears 
after the propagation. As we want to construct the largest possible sequence of monoidal 
transformations, we follow the propagation phenomenon as above. 
After more monoidal transformations: 

I ' 

: U(2,s-1) 

I , 

I 

I 

([On,m n ], [6 s+1 ,m s+1 ], [6 S , 1], [1,0], . . . , [1, 0], oo, . . . , oo) 

I 

([0„, m n ] , . . . , [5 fl+1 , m s+ i] , [9 S , 0] , [1, i] , [1, 0] , . . . , [1, 0]) 
Then, for 1 < i < s - 1, 

p(i,s) = l+p(l,s - 1) + 1 +p(2,s — 1) H h 1 +p(M - 1) 

with s — 1), 1 < I < i, defined by the induction hypothesis. 

□ 

Remark 4.6. Computation of examples in Singular with desing package has been useful to state this 
behaviour of the exceptional divisors after monoidal transformations. The implementation of this 
package is based on the results appearing in [3]. 

Theorem 4.7. Let (W,(J, c),0) be a basic object where J is a monomial ideal as in equation 
with at > c for all 1 < i < n. Then, the resolution function corresponding to (W, (J, c),0) drops 
after monoidal transformations in the following form: 
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([|,0], [1,0] 

I Xi 1 st monoidal transformation 

a 4 ? 1 , i], [i,o],..., [1,0]) 

>p (l,n) monoidal transformations 

ix t J 
([^i,0], [1,0],..., [1,0], oo) 

J. center defined only by variables in I 

([izaza,2],[i,o],...,[i,o]) 

1 1 

■ >p (2,ro) monoidal transformations 

i J 

([±^z2l, ], [1,0],..., [1,0], oo,oo) 

I center defined only by variables in I 



I 

([^,n-l],[l,0],...,[l,0]) 
I 

I 

([^,0],oo,...,oo) 



p in — l,n) monoidal transformations 



At this stage, ai > c by hypothesis, so the center of the next monoidal transformation is {Xi = 0}, 
and then we obtain an exceptional monomial. 

Proof. It follows by the propagation lemma and the fact that each time that 9 n drops E' n = Y' + 
\E\ V ^ 0, and E[ = (Y' + \E\ V ) - {E' n + ■ ■ ■ + E[ +1 ) = for all n - 1 >/> 1. □ 

Remark 4.8. Following the propagation in the previous way provides the largest branch in the 
resolution tree, because in other case, for example after the first monoidal transformation 

([^,1], [1,0],..., [1,0]) 

0], [1, 1], [1, 0] ... , [1, 0]) ([^^,2], [1,0] [1, 0]) 

looking to some chart j with j =/= i wc obtain an invariant which will appear later in the resolution 
process, after the propagation p(l, n). 

Corollary 4.9. Let (W, ( J, c),0) be a basic object where J is a monomial ideal as in equation |3J) 
with ai > c for all 1 < i < n. Therefore the number of monoidal transformations needed to transform 
J into an exceptional monomial is at most 

n-l 

1 + p(l, ri) + 1 + p(2, ri) + . . . + 1 + p(n - 1, ri) + 1 = n + ^ p(j, n). (9) 

i=i 

Remark 4.10. In this case we always have 9 n > c, so Sing(J, c) ^ at every stage of the resolution 
process. Therefore, in the resolution tree, the branch of theorem 14.71 effectively appears, and it is 
the largest, hence ([9]) is exactly the number of monoidal transformations to obtain J' = M' . 
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Proposition 4.11. Let (W, (J, c),0) be a basic object where J is a monomial ideal as in equation 
(0) with a,; > c for all 1 < i < n. Then the previous sum of propagations is a partial sum of Catalan 
numbers. 



n + n) — ^""^ Cj where Cj = ^ , — ^ ^ ^ ^ j- are Catalan numbers. 

Proof. (1) Extend p to arbitrary dimension: 

n-1 

n + ^2p(j,n) =p(n,n + 1). 

Because of the form of the recurrence equation defining p(i,j), and the fact that p{n, n) = Q 
by definition, it follows that 

n n— 1 

p(n,n + l) =n + ^2p(j,n) =n+^p{j,n). (10) 

3=1 3=1 

(2) Solve the recurrence equation defining p(i,j): 

(a) We transform the recurrence equation ([H]), defining p(i,j) for < i < j and 1 < j < n, 
to another recurrence equation defined for every j,j > 0: 

By sending the pair (i,j) to the pair (i, j — i) we extend the recurrence to i,j > 0, that 
is, we consider 

P{iJ) =p(i,i + j) 

then p(i,j) = p(i,j — i). As p(i,j) is defined for < i < j then p(i,j) is defined for 
< i < i + j for every i, j > 0. 

(b) Note that 

+ =p(i,i+j) -p(i- + 

i i — 1 

= i + J2 K fe > » + i - 1) - (* - 1) - Sp(*, * +3 ~ !) = P(*. * + 3 ~ 1) + 1 = p(*,j - 1) + 1. 
fc=i fc=i 

Therefore, we have the following recurrence equation involving p(i,j) 

P(i,j) = 1 +P( i - Ijj + 1) - 1) fori,i>l , , 

p(0,j)=p(i,0) = 1 J 

Take r(i,j) = p(i,i + j) + 1 = p(i,j) + 1 and replace p(i,j) with r(i,j) in the equation 
(jlip . It follows the auxiliary recurrence equation: 

K*j j) = r (* - 1J + 1) +r(i,j - 1) fori,j>l 
r(0, j) = p(Q, j) + 1 = 1, r(t, 0) = p(t, 0) + 1 = 1 ^ 

(c) Resolving the auxiliary recurrence equation (|12p by generating functions: 
Define r- L ,j := r(i,j) and the generating functions 

= ^ r itjX y G C[[z,y]], i2.(a;,|/) = £ r^ar'-V'" 1 G C[[x,y]]. 
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Note that R(x,y) is, by definition, the generating function of the sequence r(i,j). 
By the recurrence equation (|12[) involving r(i,j), it follows 

R s (x,y) = n-ij+ix^yi- 1 + ^ ^.j-ia^V -1 



E r^+xxV'- 1 + 5 E 

i>0j>l »>lj>0 



n,jX l y 3 



i>Oj>l 



E r «,oa; 1 + E r i,jX l y J 

i>l i>lj>l 



F E ryzV + i 



j>0,j>2 



E + E fi.^V 



E r i,jX l y J - n,ix % y 
i>a.j>i j>o 



E r o,i2/ J + xyR s (x, y) - y E ^.i 1 ' 

j>l i>0 



izj - 1 + xyR s (x,y) 

+xyR s (x,y) 



i-v 



cyR s {x,y) - y E ^,1^ 



i>0 



2/(1 



^ + br,(x, y) - \ e + + ffl»(*> y)- 



i>0 



Then 



1-y- - ) R s (x,y) 



2/(1 - y) 1 - a: y 



1 1 V- 



j>0 



multiplying the equality by y we have 



1 



(y - y 2 - x)R s (x,y) = — h 

1 - y 

+ r ,i - > r i)X x = - 

1 — 1/ 1 — x 1 



1 — y 1 — x 

y 



-y 1 - x 



i>0 



Therefore 



(y - y 2 - x)R s (x,y) 



V 



1 — y 1 — a; 



i>l 



which defines an equation of the form 

Q(x,y)R s (x,y) = K(x,y) - f/(x). 
Now apply the kernel method used in [2J, algebraic case 4.3: 

If Q(x,y) = then y = 1± v'i~ 43: . We take the solution passing through the origin, 
y = 1 ~ v/ ^~ 4a - and y = xC(x) where C(x) is the generating function of Catalan numbers. 
On the other hand, Q(x,y) — gives K(x,xC(x)) = U(x), 



K(x,y) 



y 



y 



—y 2 + y — x + 1 
1-y 1-x (l-x)(l-y) 



1 
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so K(x, xC(x)) = (1 _ I)(1 '_ IC(I)) - 1 and using jz^x) = C ( x ) we have 

Making some calculations and using that R(x, y) satisfies 

R(x,y) = xyR s (x,y) + ^r 0j V + }] r ifi x % -r ,o 

j>0 i>0 

we obtain the generating function of r(i,j) 

xyC{x) + x — y 



R{x,y) = 



(y 2 - y + x)(l - x) ' 



(3) Compute the generating function of the sequence p(n, n + 1): 
The coefficient of y in R{x, y) is just J2i>o r i,i xX then 

i dR(x,y) C{x) 

> J n,ix — 



Oy 



i=o 1 — x 



is the generating function of the elements in the first column. 

If C(x) is the generating function of C n then the convolution product C(x) ■ is the genera- 
ting function of Y^k=o ^ k ~ therefore 



r n ,i = C k . 

k=a 

As r(n, 1) = p(n, n + 1) + 1 then p(n, n + 1) = r(n, 1) — 1 = J^k=o *->k — 1, as Cq = 1 we have 

n 

p(n, n + 1) = C k 
fe=l 

where Ck are Catalan numbers. 

□ 

See [12] for more details about Catalan numbers and the web page [IT] for further details about 
their partial sums. 

Theorem 4.12. Let (W, (J, c),0) be a basic object where J is a monomial ideal as in equation {3J) 
with ai > c for all 1 < i < n. Then the number of monoidal transformations required to resolve 
(W, (J, c), 0) is at most 

n 

Cj + (2^=i °i - l)(d - c) - c + 1 

where Cj are Catalan numbers. 

Proof. It follows by theorem 12. 11 lemma [37X1 and proposition 14. Ill □ 
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Table 1: Values of the bound 



n 




global bound 


1 


i 


1 + (d - c) - c + 1 


2 


3 


3 + 7(d-c) -c+1 


3 


8 


8 + 255(d-c) - c+1 


4 


22 


22 + 4194303(d-c) - c+1 



Example 4.13. The following table shows some values of the bound for any monomial ideal J as 
in equation ([3]) with cij > c for all 1 < i < n. 

Remark 4.14. Note that, as a consecuence of proposition 14. 11[ the number of monoidal transforma- 
tions needed to transform J into an exceptional monomial only depends on n, the dimension of the 
ambient space. 

Corollary 4.15. Let J =< Z c - 1" 1 • . . . • X% n >C k[Xi, . . . , X n , Z] be a toric ideal with a* > c for 
all 1 < i < n. Then the number of monoidal transformations needed to resolve (A]J +1 , (J, c),0) is at 
most 

n 

Cj + (2 E ?=i °* - l)(d - c) - c+ 1 
where Cj are Catalan numbers and d — Yl^—i 

5 Higher codimensional case 

In the minimal codimensional case, the way in which the invariant drops essentially depends on 
the number of accumulated exceptional divisors. Because the first components of the invariant, 
6 n ,...,6\, defined in equation @, only depend on the order of the ideals I n ,...,I\. Recall that, for 
each Ji, we use the ideal (see remark 153)) . to define the ideal 7,. 

But in this case the first components of the invariant play an important role. They can also 
depend on the order of the (exceptional) monomial part M n , see remark [3. 5 1 So they may increase 
suddenly when some 6j is given by the order of the ideal M n . We will call this situation the higher 
codimensional case in dimension j. 

Note that after some monoidal transformations, we can obtain a new higher codimensional case 
in another dimension. 

So, we must compute the number of monoidal transformations while 8 n > c, with a suitable 
sum of propagations. Then, estimate the number of monoidal transformations needed to get the 
higher codimensional case in dimension 1, and use the known estimation for the order of M n to 
give an upper bound for the number of monoidal transformations needed to get the following higher 
codimensional case inside this one (if it is possible). Afterward, estimate the number of monoidal 
transformations needed to get the higher codimensional case in dimension 2, and so on. 

Hence, it has not been possible to obtain a bound for this case in the same way as above, due to 
the complications of the combinatorial problem, that perform that we can not know what branch is 
the largest in the resolution tree (to obtain an exceptional monomial). 

Furthermore, if we could find such bound, the large number of potential cases we expect, suggests 
that this bound would be very huge, even to estimate only the number of monoidal transformations 
needed to obtain an exceptional monomial. 
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Example 5.1. If we consider the basic object (W, (J, c),0) = (A|, (Xf X^X^, 4), 0), there exists a 
branch of height 15 to obtain J' = M' or Sing{J' , c) = 0. So, in dimension 3, we need a bound 
greater than or equal to 15 for a higher codimensional case, in front of the 8 monoidal transformations 
needed for a minimal codimensional case. 

Remark 5.2. In any case, both theorem l2.1l and lemma |3"TT1 are valid also in the higher codimensional 
case. So the open problem is to find a bound C to obtain an exceptional monomial, to construct a 
global bound of the form 

C +(2 C -l)(d-c) -c + 1. 

Remark 5.3. For n = 2 the higher codimensional case appears only in dimension 1 and making some 
calculations we obtain C = 3, that gives the same bound as in the minimal codimensional case. This 
bound can be improved by studying the different branches. 
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